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' Abstract 

• We present a first Wilson twisted mass fermion calculation of the matrix 

element between pion states of the twist-2 operator, which is related to the 
the lowest moment (x) of the valence quark parton distribution function in a 
pion. Using Wilson twisted mass fermions in the quenched approximation we 
demonstrate that (x) can be computed at small pseudoscalar meson masses 
down to values of order 250 MeV. We investigate the scaling behaviour of this 
physically important quantity by applying two definitions of the critical mass 
and observe a scaling compatible with the expected O(a^) behaviour in both 
cases. A combined continuum extrapolation allows to obtain reliable results 
for {x) at very small pseudoscalar meson masses, which previously could not 
be explored by lattice QCD simulations. 
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1 Introduction 



In lattice QCD the so-called chiral extrapolation of data obtained at pseudoscalar 
meson masses of about 500 MeV to the physical pion mass value (140 MeV) is one 
of the main systematic uncertainties in today's lattice QCD calculations. As an ex- 
ample from our own work let us quote a recent paper ^ where the lowest moment 
(x) of the valence quark parton distribution function in a pion has been calculated, 
fully controlling the continuum limit, finite size effects j2] and the non-perturbative 
renormalization The only remaining uncertainty, besides the quenched approx- 
imation, was the chiral extrapolation. 

In this letter we want to report on a first step towards eliminating also this sys- 
tematic error by employing the twisted mass formulation of lattice QCD |3] which 
allows to regulate unphysically small eigenvalues of the Wilson-Dirac operator and 
simultaneously achieves an 0(a)-improvement of physical observables without the 
need of improvement coefficients [Hj (automatic 0(a) improvement). With the ex- 
ample of (x) we will demonstrate that indeed this formulation of lattice QCD allows 
to bridge the gap between results at pseudoscalar meson mass values of 500 MeV 
or larger - as obtained in conventional simulations - and the physical value of the 
pion mass. Performing the simulations at a number of values of the gauge coupling 
Qq = 6//3, the continuum limit of (x) is performed at values of the pseudoscalar 
meson mass small enough to allow, at least in principle, the comparison of the 
numerical results with the predictions from chiral perturbation theory. 

Wilson twisted mass fermions have been employed already in a number of quenched 
simulations [g[7llHlinilinilIIlliaiiaiIll- See refs. O CHI for recent reviews. 

Also full QCD simulations using this approach have been performed and proved 
to be very useful in studying the phase structure of lattice QCD with Wilson type 
fermions ^3 ^1 221 120] • On the theoretical side, various studies were performed by 
means of the Symanzik expansion and of chiral perturbation theory fFR 122 12B1 1211 

l23l2nil2Z|. 

In recent quenched simulations fOl El Ej it has been shown that the Wilson 
twisted mass approach can be used to simulate small pseudoscalar meson masses of 
order 250 MeV while keeping O(a^) cut-off effects under control, when employing the 
definition of the critical mass derived from the vanishing of the PCAC quark mass 

|2ni l2Zj • These simulations were done for basic observables as extracted from 2- 
point correlation functions. In this letter we will extend the investigation of Wilson 
twisted mass fermions also to the physically important case of 3-point functions, in 
particular for matrix elements related to moments of parton distribution functions 
which are relevant in deep inelastic scattering. 
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2 Lattice action and operators 



2.1 Wilson twisted mass fermions 

In this letter we will work on a lattice xT with Wilson twisted mass fermions j3] 
which can be arranged to be 0(a) improved without employing specific improvement 
terms 0. The action for a degenerate fiavour doublet of twisted mass fermions can 
be written as 

S[U,ij,ij] = a'^^ij{x){Dw + mo + ifi'j5T3)4){x) , (1) 
where the Wilson-Dirac operator is given by 

Dw = Y. 2^^^(K + - «V;V,] (2) 
,1=1 

and V* denote the usual forward and backward derivatives, mo and fi denote 
the untwisted and twisted bare quark masses. We refer to refs. jTl |H] for further 
unexplained notations. Here and in the following ip{x) indicates a fiavour doublet 
of quarks. 

A key element in this twisted mass setup is the definition of the critical quark 
mass mc, since, in order to obtain automatic 0(a) improvement, the target contin- 
uum theory should have a vanishing untwisted quark mass. In the present work, we 
will employ two definitions for the critical mass with the final aim of a combined 
continuum extrapolation of the results obtained in the two cases. The first definition 
of the critical mass is the point where the pseudoscalar meson mass, computed with 
plain Wilson fermions (/i = 0) vanishes, the second, where the PCAC quark mass, 
computed in the Wilson twisted mass setup, vanishes. See refs. [TTJ [TTJ HH], for 
details on how the critical mass has been computed numerically. In the following 
we will refer to the first situation as the "pion definition" and to the second situ- 
ation as the "PCAC definition" of the critical point. Both definitions should lead 
to 0(a)-improvement, but they can induce very different O(a^) effects, in particular 
at small pseudoscalar meson masses [23 IIHI I2ZI • Indeed, in refs. |71 |S] we reported 
that the pion definition can have substantial O(a^) effects which are amplified when 
the quark mass becomes small and violates the inequality fi > ah? (where /i is, at 
full twist, the parameter which provides mass to the pseudoscalar meson). On the 
other hand, when the PCAC definition of the critical mass is used, these particular 
kind of O(a^) cut-off effects are dramatically reduced as was shown numerically in 
refs. [inilIIllI2], and theoretically in refs. |^ 123 123 • 

The purpose of the present letter is a demonstration that Wilson twisted mass 
fermions are in a position to reach small quark masses and eventually allow a com- 
parison with chiral perturbation theory also for more complicated physical observ- 
ables than the 2-point functions considered in refs. |nilZllHllIDllIIlinilI31IIll- A 
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computation of pion form factors employing the twisted mass fermion approach was 
aheady presented in ref. [2], but at a rather high values of the pseudoscalar meson 
mass (470 and 660 MeV), using only the pion definition of the critical quark mass. 



2.2 The twist-2, non-singlet operator 

The towers of twist-2 operators related to the unpolarized structure functions have 
the following expressions 

0;,...,A^) = ^^(a:)7{m B;., ■ ■ ■ B^,} ^r>(x) , (3) 
where {■ ■ ■ } means symmetrization on the Lorentz indices and 

D^.=D^. - D,, = ^[Vm + V;] . (4) 

The flavour structure is specified by the Pauli matrices r*^ where we include here 
also the identity with r° = 2 ■ 1. In general one should perform an axial rotation 
in order to obtain the expressions for the twist-2 operators for the twisted mass 
formulation. For our purposes it is enough to notice that the operators in eq. Q 
with flavour index a = 0, 3 do not rotate. We concentrate in this work on the 
twist-2 quark operator related to the lowest moment of the valence quark parton 
distribution function in a pion. In particular for the up quark (the down quark can 
be treated in the same way) this amounts to consider operators of the following form 

1 - (1 + r^) 
^lui.^) = 2^(^)^{m Du} ^ i^ix) - V ■ trace terms , (5) 

There are two representations of such a non-singlet operator on the lattice [2H1 
1^ . In the following we will concentrate on the operator 

Ol,{x) = ^ijix) [74 B4 D, ] ^^^^(x) , (6) 

k=l 

since in computing the matrix elements of this operator one has not to supply an 
external momentum (an external momentum increases considerably the noise to 
signal ratio). 

The matrix elements of this operator can be computed in the standard way, 
described in refs. jSOlini]. We indicate with 

P±(x) = ^(x)75^^(^) ' = ^^4^ (7) 

the interpolating operator for the charged pseudoscalar meson. The ratio of the 
3-point function 

CM = a'^5^(P+(0)O44(y,Z/4)P-(x,T/2)) , (8) 
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and the 2-point function 

Cp(x4) =a3^(P+(0)p-(x,X4)) , (9) 

X 

is related to the matrix element we are interested in. The Wick contractions of the 
correlation function (jH)) contain also a disconnected piece that we neglect consistently 
with the fact that we are interested in the valence quark distribution. In particular 
if we perform a transfer matrix decomposition and define 

R{y,) = -glii^ , (10) 

Cp(T/2) ' ^ ' 

in the limit when only the fundamental state dominates (0 -C j/4 -C T/2), and the 
ratio R reaches a plateau in ?/4, we obtain 

(0,P5|C44|0,P5) = 2mp5i? . (11) 

where |0, P^*) indicates the fundamental state in the charged pseudoscalar channel. 
The relevant bare quantity is then 

mps 

In fig. n we show an example of such a plateau from which we read off the bare 
matrix element, from our second smallest pseudoscalar meson mass and smallest 
lattice spacing. 

The matrix element obtained in this way has to be renormalized by a multiplica- 
tive renormalization factor. To this end, we took the as computed in refs. 012] 
from a Schrodinger functional (SF) renormalization scheme [321 E31 IMl with 
standard Wilson action. The renormalization factor will affect the renormalized 
matrix element by 0{a) contributions. In order to check the 0(a) effects coming 
from the boundaries of the SF we have varied the improvement coefficients typical 
of the SF boundaries and we have observed no variation in the continuum limit, 
indicating that these boundary 0{a) effects are negligible. Nevertheless the contin- 
uum limit discussed in the following section is always performed with the boundary 
improvement coefficients set to their perturbative values. 

The renormalized matrix element has a well defined continuum limit and in the 
phenomenologically relevant MS scheme it is given by 

(x) ""{jj^romps) = hm /i = 2 GeV, (13) 



where Zrgi and /'^^(/i) were computed in ref. (see this reference for further 
details). We remind here that this n = 2 GeV indicates the renormalization scale 
and not the twisted mass, since we are using a mass independent renormalization 
scale. 
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Figure 1: (x)^^'^'^ at /3 = 6.2 and a pseudoscalar meson mass of about 370 MeV; we 
show the average of the two plateaux around T/4 and 3T/4. 



3 Numerical results 

Our quenched simulations were performed for a number of bare quark masses in a 
corresponding pseudoscalar meson mass range of 270 MeV < mps < 1.2 GeV using 
the Wilson plaquette gauge action, employing periodic boundary conditions for all 
fields. In table ^ we give further details of our simulation parameters. 

We performed simulations at different values of the twisted mass parameter fi 
while setting mo to its critical value as obtained from the pion or the PCAC defini- 
tion. The corresponding critical hopping parameters can be found in ref. ^21- The 
results are summarized in table El for the pseudoscalar mass and in table El for the 
bare matrix element. 

The goal of this letter is to perform the continuum extrapolation of (x) at a fixed 
value of mpsro, for a number of values of mpsro. An interpolation of the values of 
(x) to the chosen values of mpsro is needed. These values are close to the simulated 
ones, and so even a linear interpolation is usually sufficient. By using the value of 
the force parameter rg = 0.5 fm jSHl EI] to set the scale, the lowest pseudoscalar 
meson mass that can be reached corresponds to mps = 272 MeV (for which we 
have only data obtained with the PCAC definition of Hc) ■ On the basis of the study 
performed in ref. we expect finite size effects (FSE) in the matrix element to 
be relevant for the smallest four quark masses simulated. Extending the study of 
ref. down to values of rripsL ~ 2.7 (for which the FSE can be as large as 13%), 
we have corrected the matrix elements for these effects down to the second smallest 
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/? 


5.85 


6.00 


6.10 


6.20 


6.45 


a (fm) 


0.123 


0.093 


0.079 


0.068 


0.048 


ro/a 


4.067 


5.368 


6.324 


7.360 


10.458 


L/a 


16 


16 


20 


24 


32 


T/a 


32 


32 


40 


48 


64 






pion definition 


( pion\ 

\ C J 






255 


388 


300 


207 


214 




0.0100 


0.0076 


0.0064 


0.0055 


0.0039 


Haa 


0.0200 


0.0151 


0.0128 


0.0111 




/i4a 


0.0400 


0.0302 


0.0257 


0.0221 




H^a 


0.0600 


0.0454 


0.0385 


0.0332 




Hea 


0.0800 


0.0605 


0.0514 


0.0442 




jija 


0.1000 


0.0756 


0.0642 


0.0553 








PCAC definition 


C /.PCAC ) 
\ c J 




N 


400 


300 




300 




jjLia 


0.0050 


0.0038 




0.0028 




H2a 


0.0100 


0.0076 




0.0055 




fisa 


0.0200 


0.0151 




0.0111 




jjL^a 


0.0400 


0.0302 




0.0221 




Hr^a 


0.0600 


0.0454 




0.0332 




/iga 


0.0800 


0.0605 




0.0442 




jjLja 


0.1000 


0.0756 




0.0553 





Table 1: Simulation parameters and the number of measurements (N^cas) 

quark mass (mps=368 MeV). For the smallest quark mass, however, the sensitivity 
required to investigate FSE is computationally very expensive and here we present 
the corresponding points without corrections, with the purpose of showing that, even 
for quantities more complicated to extract than meson masses or decay constants, 
there are no problems of principle in reaching small quark masses. 

In fig. 121 we show the combined continuum extrapolation of (x) obtained with 
the two definitions of already converted to the MS scheme at /i = 2 GeV as 
explained in the previous section, for a wide range of values of fixed pseudoscalar 
meson masses. In principle, employing the renormalization factors obtained with 
(untwisted) Wilson fermions 0(a) lattice artifacts can be introduced, which are ab- 
sent in the bare matrix elements. However this kind of 0(a) effects are independent 
of the mass and of the definition of the used. Considering for example the case 
of the second lowest mass /i2 {rrips = 368 GeV) in fig. |21and performing a combined 
continuum fit of the type A + Ba/vQ + C(a/ro)^ where A and B are the same for the 
two definitions of Hc while C is different, one finds that B ^ 0. Since B/A does not 
depend upon the mass it follows that in the determination of the renormalization 
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5.85 


6.00 


6.10 


6.20 


6.45 








0.2240(33) 


0.1773(43) 


0.1482(27) 


0.1282(24) 


0.0892(22) 


fi3a 


0.3117(23) 


0.2379(29) 


0.2030(21) 


0.1760(20) 




jji/^a 


0.4430(33) 


0.3337(22) 


0.2865(15) 


0.2472(16) 




H^a 


0.5523(23) 


0.4135(17) 


0.3534(13) 


0.3055(14) 






0.6478(21) 


0.4840(16) 


0.4130(13) 


0.3561(12) 




fija 


0.7349(21) 


0.5491(14) 


0.4676(12) 


0.4034(11) 








Hia 


0.1640(24) 


0.1178(69) 




0.0936(22) 






0.2289(19) 


0.1686(50) 




0.1278(20) 




/isa 


0.3231(14) 


0.2401(33) 




0.1781(17) 




^40 


0.4608(12) 


0.3422(23) 




0.2495(13) 




li^a 


0.5703(11) 


0.4233(18) 




0.3080(11) 




l^ea 


0.6659(10) 


0.4943(15) 




0.3588(10) 




jjL-ja 


0.7532(9) 


0.5595(15) 




0.4062(9) 





Table 2: Pseudoscalar meson masses mpsa for all simulation points. These data 
refer to a subset of the data obtained in ref. liif . 

factors 0(a) lattice artifacts are in practice negligible. We are thus justified in per- 
forming a continuum extrapolation of the type A + C(a/ro)^ and we can see from 
fig. I21that the scaling of (x) is in agreement with pure O(a^) cut-off effects for both 
definitions of the critical mass and values of /3 > 6.0. The slope of (x) as a function 
of appears to be rather small for the PCAC definition of Kc- At /3 = 5.85 and 
lower, lattice artifacts which increase with /i are visible for the highest masses while, 
in analogy to what observed in ref. jT2], by using the pion definition of Kc O(a^) 
cut-off effects are enhanced at small quark mass. Indeed, for the second smallest 
pseudoscalar meson mass (the smallest one is absent with the pion definition of Hc), 
we performed an additional simulation at /5 = 6.45, in order to have a better control 
on the continuum extrapolation, and we excluded the point at /? = 6.0 (which ap- 
pears to be outside of the scaling region). Fig. |21 nicely demonstrates that using only 
this definition of the critical mass it is important to add the data point at /5 = 6.45 
and this effect would be probably even worse at the smallest pseudoscalar meson 
mass. 

In fig.|21and table|3]we present the results for {x)^^{fi = 2GeV) in the continuum 
as a function of the pseudoscalar meson mass in GeV. The empty squares are our 
values obtained earlier from a combined continuum extrapolation of Wilson and 
clover-improved Wilson fermions data using the Schrodinger functional scheme 
As usual, such quenched simulations have to stop at a pseudoscalar meson mass of 
about 600 MeV. At such high masses it becomes very difficult, if not impossible, to 
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5.85 


6.00 


6.10 


6.20 


6.45 




0.90(5) 


0.95(4) 


0.99(4) 


1.01(4) 


1.06(5) 








0.2848(105) 


0.2143(107) 


0.2028(93) 


0.2073(103) 


0.2226(108) 


fi^a 


0.3295(57) 


0.2818(56) 


0.2700(53) 


0.2621(49) 






0.3643(37) 


0.3294(33) 


0.3144(28) 


0.2996(33) 






0.3840(25) 


0.3533(24) 


0.3374(21) 


0.3215(26) 






0.4012(20) 


0.3711(20) 


0.3551(17) 


0.3386(24) 






0.4168(18) 


0.3861(17) 


0.3700(15) 


0.3539(20) 








Hia 


0.2566(157) 


0.2615(219) 




0.2505(241) 




ii2a 


0.3049(80) 


0.2819(110) 




0.2698(98) 




H3a 


0.3420(39) 


0.3135(58) 




0.2907(46) 




jji/^a 


0.3704(24) 


0.3439(36) 




0.3127(31) 




fi^a 


0.3900(18) 


0.3631(26) 




0.3300(24) 




/ieo 


0.4077(15) 


0.3791(21) 




0.3455(19) 






0.4237(13) 


0.3935(18) 




0.3502(16) 





Table 3: Renormalization factor = ZRGi(a)//^^(/i) for fi = 2 GeV from ref. |7]/ 
and bare matrix element (x)^^^'^ for all simulation points. 



mps [GeV] 




0.272* 


0.260(31)* 


0.368 


0.243(21) 


0.514 


0.272(21) 


0.728 


0.299(22) 


0.900 


0.317(23) 


1.051 


0.335(24) 


1.163 


0.350(25) 



Table 4: (x)^^ in the continuum using a combined extrapolation of data obtained 
with the PCAC and pion definition of k^* The value corresponding to a pseudoscalar 
meson mass of 272 Me V has not been corrected for FSE. 

compare the simulation results to chiral perturbation theory [3H1 123 HO] or to other 
phenomenological predictions |H1 132] , even when the results are extrapolated to the 
continuum limit as done here. 

Fig. 121 shows that with Wilson twisted mass fermions, the large gap between 
pseudoscalar meson masses of about 600 MeV, as the lower bound for standard 
simulations, and the physical value can be bridged. Quenched chiral perturbation 



8 



theory predicts the absence of chiral logs for the matrix element studied in this 
letter, but our present large error bars, and the lack of more data in the region 
('^ps ^ 500 MeV) where chiral perturbation theory should be applicable, does not 
allow us to perform a careful chiral extrapolation. We quote then as our final result 

(a;)^(2GeV) = 0.243(21) (14) 
given by the value at the next to smallest pion mass. 

4 Conclusions 

In this letter we performed a test of Wilson twisted mass fermions for interesting 
physical observables, the moments of parton distribution function which are relevant 
in deep inelastic scattering. So far, most investigations of Wilson twisted mass QCD 
at small masses considered only 2-point correlation functions. The present study 
is the first that investigates 3-point correlators down to masses of order 250 MeV. 
In particular, here we have studied the example of a twist-2 operator. The matrix 
element of such a renormalized operator between pion states (x) corresponds to the 
average momentum of the valence quark distribution {up for example) in a pion. 

In the present work we employed two definitions of the critical mass, the pion 
and PCAC definition fUE]. The scaling of (x) is in agreement with the expected 
O(a^) cut-off effects for both definitions of the critical mass. We could perform a 
controlled continuum extrapolation of (x) down to pseudoscalar meson masses of 
about 270 MeV by combining the data obtained with the two definitions. Of course, 
in principle, also overlap simulations are able to reach such values of the pseudoscalar 
meson mass. This will come, however, at a much higher simulation cost Our 
final figure, fig. El clearly demonstrates that with our present setup it is possible to 
bridge the gap between large pseudoscalar meson mass values of 600 MeV and the 
physical value of the pion mass. 
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Figure 2: {x)^^ as a function of {a/ro)^ for different values of the pseudoscalar 
meson mass. 



14 



0.5 
0.45 

0.4 - 
0.35 - 

0.3 
0.25 - 

0. 
0.15 

0.1 - 
0.05 - 




' 1 

<X>'^^ (^i=2 GeV) 



0.2 -{ 



• Wilson tmQCD at nil 
□ combined Wilson-Clover 
<^ experimental value 



mj [GeV'jJ 



0.8 



1.2 



1.6 



Figure 3: {x)^^{fi = 2 GeV) extrapolated to the continuum as a function of the 
pseudoscalar meson mass. Squares are obtained from a combined continuum extrap- 
olation of earlier Wilson and clover improved Wilson simulations p. The circles 
represent our results using Wilson twisted mass fermions. For the empty circle at 
the smallest mass see the text. The diamond represents the experimental value as 
obtained from global fits |^ E3] . Recently a new analysis [IHj has been performed 
giving as a result = 2.28 GeV) = 0.217(11). 

°Wc thank C. Roberts to bring this reference to our attention. 



15 



